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Abstract 



We consider application of time-delayed feedback with infinite recursion for control of anharmonic (nonlinear) oscillators subject 
(— | to noise. In contrast to the case of a single delay feedback, recursive delay feedback exhibits resonances between feedback and 
O nonlinear harmonics, leading to a resonantly strong or weak oscillation coherence even for a small anharmonicity. Remarkably, 
& these small-anharmonicity induced resonances can be stronger than the harmonic ones. Analytical results are confirmed numerically 
C for van der Pol and van der Pol-Duffing oscillators. 
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1. Introduction 

Delayed feedback is an effective control tool for chaotic 
and stochastic dynamics of solely oscillatory systems yj, y, 
and collective dynamics of ensembles of os- 
cillators S H [HI ll 111 . For a recursive delay feedback (or 



"extended delay feedback") the forcing on the system at the 
time moment t is determined not merely by the state difference 
z(f) = [x(f — t) - x(f)] (t is the delay time) but by a chain 
of differences z(f) = 2+!° fl"[x(f - (« + 1)t) - x(f - tit)} with 
\R\ < 1 Il2lll3l1 . The employment of such a feedback can pro- 
vide new opportunities in comparison to the case of the "single 



delay feedback" control IU2L Il3l Il4i 11511 . Moreover, the feed 



back has such a form when the system is connected to a res- 
onator (for instance, an optical delayed feedback for lasers is 
realized with a Fabry-Perot interferometer). A recursive delay 
feedback is also constructively inherent to some electronic de- 
vises (e.g. HH]). It is noteworthy that, from the mathematical 
perspective, the infinite recursion is equivalent to the addition of 
the neutral equation s(f) = Rs(t - r) + x(f - t) (then the feedback 
signal z(f) = s(f) - x(f))- Properties of this neutral equation 
provide favourable conditions for resonant phenomena in the 
system. 

For limit-cycle oscillators subject to noise, the regularity of 
dynamics is characterized by oscillation coherence [17] and 
"reliability" 111811 . Coherence can be quantified by the diffu- 



sion constant of the oscillation phase ip(t). While for van- 
ishing noise the phase grows linearly with time, noise results 
in random excursions of the phase around the average growth 
trend, ([<p(t) - ip(0) - (<p)t] 2 ) « Dt, where D is the diffusion 
constant. The "reliability" is the ability of a system to pro- 
vide one and the same response for the noise signal of a pre- 
recorded waveform. Mathematically, reliability means stability 
of a system response and can be measured by the Lyapunov 
exponent lfl9l 2oL 21, 22, 23]. The effect of a control on the 
dynamics regularity can be, therefore, tracked with the phase 
diffusion constant and the Lyapunov exponent. 

The employment of single and/or recursive delay feedback 
control received previously analytical treatments only for the 
case of a harmonic oscillator |^|7lll5|]. In this Letter, the analyt- 
ical treatment is generalized for anharmonic oscillators. Strong 
anharmonic resonances are revealed even for a relatively small 
anharmonicity. 

The paper is organized as follows. In Sec. |2]we derive alge- 
braic equations for the mean frequency, the diffusion constant, 
and the Lyapunov exponent for general limit-cycle oscillator 
subject to noise and a general linear feedback represented in 
terms of a Green's function. In Sec. |3]we consider the form of 
these equations specific to a recursive delay feedback and treat 
the role of anharmonicity. The analytical theory we construct is 
applied to van der Pol and van der Pol-Duffing oscillators. The 
anharmonic resonances are also observed with a direct numer- 
ical simulation. Finally, conclusions are summarized in Sec. [4] 
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2. Phase description 

Let us consider an A^-dimensional limit-cycle oscillator sub- 
ject to recursive delay feedback and noise: 



*; = Fi(x) + a Zi(t) + Bj(x) o 



(1) 



where i = 1,2, ...,N, a is the feedback strength, azi(t) is the 
feedback term, £(f) is white Gaussian noise, (£) = and 
<f(0£(O> = 25{t - f), "°" indicates the Stratonovich form of 
stochastic equation. For briefness, we will first keep the feed- 
back term in a general form valid for any linear feedback (as in 
Refs. HIH) 

+ " N 

Zi(t)= \ YjGijitOxjit-t^dh, 
o ■/=' 

involving Green's function Gjj(t); for recursive time-delayed 
feedback 

G(f) = 2S JV[<5(* - (« + 1)t) - 5(t -m- 0)] 

n=0 



= 2S 



- S(t -0) + (l-R) YR"- [ 6(t - tit) 



n=\ 



(2) 



where S is a constant matrix, t is the delay time, \R\ < 1, and 
we explicitly indicate that JT 6(t - 0) dt — 1 . 

For a noise-free limit-cycle oscillator, the oscillation phase 
<p = (/?(x) can be introduced on the limit cycle xo(f) = xo(f + To) 
and within its finite vicinity in such a way that ip — Qo = 2n/To, 
where Q.o is the natural frequency of the oscillator 11261 12711 . The 
limit cycle can be parameterized with the phase, XoQp) = Xo(<P + 
2n). In the presence of a weak forcing (noise and feedback) the 
phase description can be still utilized; 



N N 



= Q. + a J ^^Gij(ti)Hij(<p(t-ti),<p(t))dti 

+e/(v(0)°f(t), 



'=1 7=1 



(3) 



where 



j = l \ UX J tXrXobpi) 



is a 27r-periodic function featuring the sensitivity of the phase 
to noise, s is the noise amplitude, //, ; (i/f, <p) is the increase of 
the phase growth rate created by the feedback term Xj(if/) acting 
on the variable Xiifp); 



Hn(ifr,<pi) :-- 



d<p(x) 



dxi 



X=Xo((Pl) 



Notice, accurate derivation of Eq. © ||28l |29[ |30| yields a drift 
term oc s 2 , which represents the role of the amplitude degrees 
of freedom. Since Eq. (0 provides a noise-induced mean fre- 
quency shift oc e 2 , i.e., of the same order of magnitude, this 
omitted term should be accounted for, when one considers the 



effect of noise on the mean frequency. However, for robustness 
quantifiers (the phase diffusion constant and the Lyapunov ex- 
ponent, which are of our interest here) this term has been shown 
to be negligible POT . 

For weak noise, the linear in noise approximation is relevant 
and yields (see Appendix for detail) the mean frequency 



Q = C2 + a I Yj^jGijiOhiji-nOdt, 



(4) 



i=l 7=1 

the phase diffusion constant 

2s 2 (f% 



D 



N N 2 

(l+a J t 2 2 G u (t)h' {-Q.t)dt) 

'=1 7=1 



= 2s 2 (f 2 ), 

and the leading Lyapunov exponent 

s 2 {{f')\ 



3D. 

xTo 



(5) 



+°° N N 2 

(l+a / t 2 Z G u (t)h'(-£lt)dt) 

1=1 7=1 



s 2 {{f)\ 



(6) 



where (...^ = {In) 1 J Q " ...dtp, 

hijiifi) := (flijifp + tff, 

is the average susceptibility of the phase to the feedback term 
Xj(<p + i//) acting on Xj(tf) (notice, this means a phase delay by 
-i/f), and the prime denotes derivative. Eq. (0) can be solved 
with respect to O, and then, with O evaluated, Eqs. (0 and (|6]) 
provide D and /I. Noticeable constancy of the ratio 

-A = <(/0 2 ) 

O 2(/ 2 ) 

was aheady discussed in the literature ||7l|25l0]. Henceforth, 
we do not consider A because considering the diffusion constant 
is enough. 



3. Recursive delay feedback control for harmonic and an- 
harmonic oscillators 

3.1. Harmonic oscillator 

Previous analytical considerations were focused on the case 
of a nearly harmonic oscillators |H1 HI, 24], such as the van 
der Pol one; 



x = y, y = fi(l - 4x )y - x + az(t) + eij(t) 



(7) 



with /j <<c 1 (/i characterizes closeness to the Hopf bifurcation 
point). For/i <K 1, the phase <p = -arctan(y/x) and the limit 
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Figure 1; Van der Pol oscillator l|7] with /i = 0.7 subject to the recursive delay 
feedback control, (a): The phase portrait of the system without noise and con- 
trol, (b): The function of susceptibility to feedback ^22(<A) (black curve) and its 
harmonic part (blue curve), (c): The analytical dependence of the phase diffu- 
sion constant on the delay time r (black curve) and the harmonic-approximation 
version of this dependence (red curve) for R = 0.5 and a = 0. 1 (Do and 
Tq ss 2^/0.97 are the diffusion constant and the mean oscillation period of the 
control-free system, respectively), (d): Results of a direct numerical simulation 
with R = 0.5, a = 0. 1, noise strength s = 0.02. 



cycle is: xo = cos ip, yo = - simp; therefore, the equations for 
H u {\p, ip) and f(ip) yield 



Figure 2: Van der Pol oscillator l|7] with fi = 1,0 subject to the recursive delay 
feedback control, (a): The phase portrait of the control-free noiseless system, 
(b): /i 2 2(iA). (c): Analytical dependence D{t) (T ss 2^/0.94). (d): Numerical 
simulation. For detail see caption to Fig.^ 



The recursive delay feedback az(t) is typically implemented in 
the form 



z(t) = 2 2 R"(y(t - (n + 1)t) - y(t - nr)) , 



(8) 



which means 



h(<A) = 



■ sin iff j cos if/ 
i cos iff A sin (// 




1 



f(<p) = -sin<p. 



as we have only the y- variable acting on y (compare Eq. O and 
Eqs. (Q}, ©), and with hn(^i) = (1/2) sin ifr one obtains for a 
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Figure 3: Van der Pol-Duffing oscillator 114) with /i = 0.7 and b = 0.5 subject 
to the recursive delay feedback control, (a): The phase portrait of the control- 
free noiseless system, (b): /i220W- ( c ) : Analytical dependence D(r) (To ~ 
2ff/1.14). For detail see caption to Fig.^ 



Figure 4: Van der Pol-Duffing oscillator 1141 with /j = 1.0 and b = 0.5 subject 
to the recursive delay feedback control, (a): The phase portrait of the control- 
free noiseless system, (b): hiity). (c): Analytical dependence D(t) (Tq k 
2^/1.12). For detail see caption to Fig.^ 



harmonic oscillator (cf 1 15]) 



(1 -fl)sinQr 



D = 



1 + R 2 -2RcosQt ' 
D 



1 +ar(l -R) 



(1+R 2 )cosQ.t-2R 
(1 +fl 2 -2flcosQr) 2 



(9) 
(10) 



where Do = 2s 2 (f 2 )^ is the phase diffusion constant of the 
control-free system. 

3.2. Anharmonic oscillator 

Generally, oscillators are anharmonic and /i22(<A) has addi- 
tional terms, not only sin xjj. Let us consider /i22(<A) in the form 
of a Fourier series; 

oo 

hnW = £j{a m sinm^ +J3 m cos mi//) . (11) 

m— 1 

Notice, the term [3\ cos iff belongs to the first harmonic of the 
series as well as a\ sini/r, but it is owed to anharmonicity; for 
harmonic oscillations /i22(<A) is purely proportional to sin tfr. 
With Eq. (Qj}, Eqs. © and © take the form 



Q = Q 

OO 

m= 1 



a m (l - R) sinrnQr + /3 m (l + R)(l - cos wOt) 
1 + R 2 - 2R cos mQ.T 



, (12) 



D = Di 



1 +2arar(l - R) 



m=\ 



a m [(l + R 2 ) cos rnQi - 2R] + /3JI - R 2 ) sin wOr 



(1 +R 2 -2R cos mQr) 2 



-2 



(13) 



3.3. Anharmonic resonances 

Denominators in the sums in Eqs. (TT2l and ([TBI have extreme 
values for cos mQ.T — +1. For a single delay, R — 0, they equals 
1, while even for moderate values of R the minima of the de- 
nominator can be close to zero. Remarkably, at points of these 
resonances (cosmQr = +1) contribution of the mth harmonic 
can be large even when a m and /3 m are relatively small. For 
instance, next to t such that |cos3£2r| = 1, but |cos£2t| + 1, 
one can expect resonantly strong or weak phase diffusion, due 
to anharmonicity. 

For the demonstration and invigoration of the analytical find- 
ings, we have considered a van der Pol-Duffing oscillator: 

x = y, y = yu(l - 4x 2 )y - x - bx 3 + az(t) + sg(t) , (14) 

where b is the Duffing parameter, b — corresponds to a van 
der Pol oscillator ((7j ; feedback term az(t) is given by Eq. ©. 

For this oscillator, hnW) and its Fourier transform dm were 
calculated. Figs. [T]© present samples of limit cycles (a) and 
functions /122 W (b); the corresponding Fourier coefficients a m 
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b = 0, ij = 0.7 


b = 0, n=1.0 


b = 0.5, n = 0.7 


6 = 0.5, //=1.0 


a,n 


Pm 


a m 


Pm 




Pm 


a,n 




m — 1 


+0.501700 


-0.076145 


+0.499449 


-0.093296 


+0.506334 


+0.153265 


+0.505712 


+0.070064 


m — 2 


























m — 3 


-0.010249 


-0.003886 


-0.018180 


-0.010697 


-0.009246 


-0.001205 


-0.015938 


-0.005048 


m — 4 


























m — 5 


+0.000044 


+0.000084 


+0.000035 


+0.000422 


+0.000050 


-0.000003 


+0.000149 


+0.000107 



and p m are provided in Tab. [TJ (Interested readers can find a 
Maple program with comprehensive calculations in the supple- 
mentary material.) One can see in the table, that typical values 
of higher harmonics are small: and are only a few percent 
of a\ « 1 12. Nevertheless, the dependence D(t) exhibits large 
peaks related to the m = 3 resonances (in Figs. QJI4] these de- 
pendencies are plotted for R = 0.5). The results of calculations 
without higher harmonics are plotted with red curves and can 
be compared to the one of the "all-harmonics" calculation; res- 
onant discrepancies can be clearly seen. In particular, the m — 3 
resonances are even stronger than the harmonic ones related to 
m = 1 . The results of a direct numerical simulation show these 
peaks as well (Figs. QJ and[2jl) and are in a good agreement 
with the analytical dependencies. 



4. Conclusion 



We have considered limit-cycle oscillators subject to weak 
white Gaussian noise and recursive delay feedback control. In 
this way, we have derived the phase reduction equation and 
calculated analytically the mean frequency, the phase diffusion 
constant, and the Lyapunov exponent (Eqs. ©-(O) for the case 
of a general linear feedback ("single" and recursive delay feed- 
back, linear frequency filter, etc.) and a general limit-cycle os- 
cillator. 

It has been found that even a small anharmonicity leads to 
additional resonances and resonantly strong or weak coherence 
(measured by the phase diffusion constant). For instance, for 
van der Pol and van der Pol-Duffing oscillators, the anhar- 
monic resonances due to the harmonic ft™ (i^) = aj, sin 3ijr + 
/?3 cos 3i[/ are even stronger than the harmonic resonances due 



to h^^iiff) = (l/2)sinifr (Fig. [TjQ. Hence, for application 
of the recursive delay feedback control, anharmonic resonances 
are not to be neglected even for nearly harmonic oscillators. 



Appendix A. Derivation of Eqs. ©, (|5), and © 



The derivation presented here is performed in the same way 
as in Refs. ^ for the diffusion constant and in Ref. [7] for the 
Lyapunov exponent. Eq. Q can be rewritten in terms of ip — 



(A.l) 



Q. + v(i), where Q, is the mean frequency and (v(f)) = 0; 

dt\ V Gij(h) 
Tj 

/t-ti nt 
v(t 2 )dt 2 , Qf + I v{t 2 )dt 2 ) 

+ e/(Ot+ J v(fe)<*fe)of(f) 

The mean frequency and robustness quantifiers are related to 
long-term behavior of v(f) and, therefore, we can perform av- 
eraging with respect to "fast" variable Q.t (rigorous multiscale 
analysis yields the same final expressions but in a very lengthy 
way and is not presented here). Eq. (1A.U turns into 

Q. + v(f) = Qo + a \ dt\ 2 ' Gjj{t{) 



X 



X+oo 
dh^Gtji 

i.j 

hJ-Qh - f v(t 2 )dt 2 ) 
Jt-tt 

+ sf(nt + J v(t 2 )dt 2 )°t(f), 



(A.2) 



where fty(^) := (tfyftc + ip,<p))<p, (...)<p = (2n) 1 J Q " ...dtp. Lin- 
ear in v approximation reads 

X+oo 
dt\ ^ Gij(t{) 

x[Ay(-nfi)-A;/-»i) f v(t 2 )dt 2 ] 

[f(Qt)+f'(nt) J v(t 2 )dt 2 ]o£(t). (A3) 



Averaging Eq. (IA.3b yields 



y 

Ti 



Gij(t) hij(-Qt) dt + 0(e 2 ) . (A.4) 



d. = a 



Further, multiplying Eq. (1A.31 by f(f + s) and averaging, we 
can evaluate (£(f + s) y(f)) from known (£(f + s) £(t)). Multiply- 
ing the same equation by v(f + s) and integrating, we evaluate 
(v(f + s) v(f)) from (v(f + s) £(/))■ The integral of the autocorre- 

X+OO 
(v(f + s) v(f)) ds, yields the phase diffusion 

constant D; 



D 



2s 2 (f%(l+0(s 2 )) 



(l+a j^°tZGtj(f)^-CU)dt) 



(A.5) 
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Similarly to Ref. [7] we find that the leading Lyapunov ex- 
ponent 

s 2 ((f') 2 )Al+0(s 2 )) 

A = - ^ ■ ( A - 6 ) 

(l+aj^tZGijiOh'^-nOdt) 

u 

Eqs. dA.41 >— dA.6b with s <k 1 turn into Eqs. dUi-©. Remark- 
ably, the ratio A/D is constant only in the leading order with 
respect to e; the proportionality between A and D does not hold 
for higher order corrections. 
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